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Abstract 



We perform a covariant wavepacket analysis of neutrino oscillations taking into ac- 
Q I count the lifetime of the neutrino production process . We find that flavour oscillations 

Q"^ I in space are washed out when the neutrinos are produced from long lived resonances 

^ I - and what may be observed in appearance/disappearance experiments is a uniform 

>- I conversion probability independent of distance. The lifetime of the resonance which 

produces the neutrinos acts as the the effective baseline of the experiment. For this 
"^ I reason the LSND experiment where neutrinos are produced from muon decay has two 

Q I orders of magnitude more sensitivity to neutrino mass square difference than other ex- 

^ I periments where the neutrinos are produced from pion or kaon decays. Wea also show 

■^ I that there are no EPR type oscillations of the secondary decay particles. 

^. ■ 

•-H ; 1 Introduction 

d I There are two different formulas which describe flavour oscillations in spacetime; the one applied 
to neutrino oscillations |l|] is derived in the high energy regime and the other applied to kaon 
oscillations 0] is valid in the non-relativistic limit. From the derivations of these formulas it is not 
clear what formula is applicable at intermediate energies (for example for kaon oscillations from 
stopped protons where [m? /E'^) ~ 25%) . A covariant derivation of a flavour oscillation formula 
which would be valid at both high and low energies would be of interest from the conceptual as 
well as the experimental point of view. 

Kayser and Stodolsky have advocated a covariant generalisation of the non-relativistic 
phase factor exp{—imt) by replacing the absolute time t by the Lorenz invariant proper time 
s = (t^ — x^)^/^. In the lab frame the covariant expression for the phase factor is exp — irriiSi = 
exp {—imjt/Ei). The phase difference between two mass eigenstates is A[rn? / E) ~ {Am'^/E) which 
is twice the the phase difference of the standard formula |1|. This result has prompted the claim 
1^ that kaon oscillations in $ factory will have oscillation length which is half of what is given by 
the standard formula. This claim has been refuted by ||^ who take the view that the interference 
phase difference should not be evaluated at different space-time points but should be evaluated at 
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the average spacetime interval and taking the phase difference to be (mi — m2)(si + S2)/2 instead 
of {rriiSi — 771282) they recover the standard formula. Other methods of showing that an extra fac- 
tor of two does not appear in the kaon oscillations formula have been discussed in |^. A covariant 
derivation of the neutrino oscillation formula has been given by Grimus and Stockinger ||^ who treat 
the entire process of neutrino production ,propagation and detection as a single Feynman diagram. 
They show that on taking the large distance limit of the neutrino propagator the scattering cross 
section shows a space-time oscillatory behaviour and the oscillation length is identical to that given 
by the standard formula [|ll. G-S assume the initial states to be plane waves therefore the concept 
of a coherence length |^ does not emerge in their formulation. 

The oscillation phenomenon occurs because particles are produced and detected as weak in- 
teraction eigenstates (the neutrino states u^, z/^ and z/^- or the Kaon states K^, K^ etc) but the 
propagators are diagonal in the mass eigenstates (the neutrino mass eigenstates z/j, -i = 1 — 3 or the 
Kaon mass eigenstates Kl,Ks)- The probability amplitude for oscillations of a gauge eigenstate 
(|a >) to another (|/3 >) is a lof a gauge eigenstate (|a >) to another (|/? >) is a linear superposition 
of the propagation amplitude of the mass eigenstates {\i >). 

Aia ^ (3-t)=Y, <P\i> < i\e-'"% ><i\a> (1.1) 
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In relativistic field theory the propagation amplitude of the main relativistic field theory the prop- 
agation amplitude of the mass eigenstates < i\e~^^^\i > can be identified with the Feynman prop- 
agator. At large time-like spacetime separation, the Feynman propagator (for both fermions and 
bosons) has the form, 

where Si is the invariant space-time interval propagated by the z/j mass eigenstate. In the lab 
frame the phase factor of the amphtude ( |1.2D is given by exp ( — ^^ ) .In this paper we derive the 
oscillation amplitude by evaluating the Feynman propagator of the neutrinos in the large time-like 
asymptotic limit. The asymptotic propagator of a position eigenstate has the form K{x,t]7ni) ~ 
{mi/2nisi)^^'^ exp—imiSi where Sj is the spacetime interval propagated by the 77ii mass eigenstate. 
We derive the same expression in the path integral method. Using the path integral method we show 
that for large distance propagation the forward in time paths are equivalent to on-shell particles. 

Another aspect of flavour oscillations on which there is no consensus of opinion is the ques- 
tion of the energy momentum of a flavour eigenstate. A flavour eigenstate propagates as a linear 
combination of different mass eigenstates which are on-shell. Some authors have argued that the 
different mass eigenstates have common energy but different mommenta Pi = \ E'^ — 7nf while 



others have assumed that the states have same mommenta but different energy Ei = JP^ + rnf. 
We show here that different mass eigenstates in a linear combination of states have niether energy 
nor mommentum in common. For long distance propagation, the energy and mommenta of a given 
mass eigenstate are not independent however but are related by the mass shell condition. We use 
these results to show that there are no EPR type associated particle oscillations as has been claimed 
in literature 0. This result is in agreement with the conclusions of Lowe et al. and Dolgov et 
al. 10 . Our proof however is more direct in that we show that although the energy and momenta 
of the muons depend on the outgoing neutrino mass, the phase difference of the muon states cor- 
ressponding to different neutrino masses cancels exactly and there is no spacetime dependence of 
the muon probability distribution. 
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If the initial wavef unctions of the propagating particles were strictly delta functions in spacetime 
then no interference between different mass eigenstates can take place. We therefore generalise the 
delta function propagators to propagators of gaussian wavepackets. The interference term as a 
function of distance is obtained by taking the time-overlap of different mass eigenstate propagators 
. The expression for flavour conversion probability (for say two flavours with mixing angle ^) as a 
function of distance X turns out to be, 

P{a ^P;X) = i W2^(l - cosi^X) e'^) (1.3) 

Therefore extra factor of two which came from naively subtracting the phases of plane wave prop- 
agators 1^ goes away in the wavepacket averaging and the standard expression for the oscillations 
length is recovered. This is in agreement with [§-0. The covariant wavepacket treatment intro- 
duces a new contribution to the exponential factor A when particles are produced from long lived 
resonances ( for example for neutrinos from muon decay as opposed to neutrinos from Z decay) . At 
distances smaller than the coherence length , the exponential suppression factor is, 

, Aim? ,0 , , 

where r is the lifetime of the resonance which produces the particle which undergo flavour oscilla- 
tions. If the uncertainty in position of the initial particle vt is larger than the detector distance 
X , the exponential term washes out the oscillations in ( |1.3| ). In neutrino experiments where the 
source is pions, kaons, muons or nuclear fission, the spatial oscillations of the conversion probability 
cannot be observed. What can be observed in these experiments is a constant (distance indepen- 
dent) conversion probability. The conversion probability is sensitive to values of Am^ ^ {2\/2E/t) 
, which means that the experimental bound on Am^ is lower with longer -lifetime sources. For 
this reason the LSND experiment [^ which uses neutrinos from muon decay (r^ = 2.19 x 10~^s) is 
sensitive to two orders of magnitude lower neutrino mass square difference compared to other accel- 
erator experiments like BNL-E776 |jlO| , Karmen [|lT] and CCFR |]12| which use neutrinos from pion 



and kaon decays (r ~ 10~*sec). We fit the experimental data from LSND along with BNL-E776, 
Karmen, CCFR and Bugey [|1^] experiments with the conversion probability formula ( |1.3| ) and plot 
the allowed range for the mass difference and mixing angle. We find that there is a large region 
of LSND allowed band which is not ruled out by the other experiments. The region of Am,^ and 
sin'^26 which is allowed by LSND and all other experiments is 2.2 x lO^'^el^^ < Am^ < Q.3eV^ and 
0.002 < 3171^29 < 0.15. The lower bound on allowed Am? is two orders lower than what is obtained 
by fitting the LSND data with the standard oscillation formula (equation (|1.3| ) with A set to zero). 

2 Field theoretic derivation of phase factor 

In relativistic field theory the propagation amplitude of the mass eigenstates < i|e~*^*|z > can 
be identified with the Feynman propagator Spixf — Xi,mi) =< T Vi{xi)i'i{x j) >. The Feynman 
propagator in position space is , 

K(x,mi) = -I / - — T7-^-^ — o — ^e *^''. (2.1) 

^ ' ^ J {2-kY p^ - m} - le ^ ' 
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This integration can be done by expressing the denominator as an exponential, 



-7, n = / da exp \ia(p'^ — mf — ie)} (2.2) 

p^- mi -It Jo 



and integrating over the resulting Gaussian in p. The remaining integral over a, 

/■OO ('^2 |xP) 

/ da q;~^ exp {iami + i } (2.3) 

JO 4a 

can be performed by substituting s = (t^ — |xp)^/^, z = ims and i] = 2(am/s) and making use of 
the integral formula [Il| for the Bessel function 



1 roo y 1 

- / dr] r/-(''+i) exp - -(r^ + -) = K,{z). (2.4) 

/JO Z r] 

The resulting expression for the Feynman propagator ( p.l|) is 

? Tn ■ 

K{x,m,) = {-—){t ^ + m,) (^) K,{tm,s,) (2.5) 

47r^ Si 

where Sj = (t^ — x^)^/^ is the invariant spacetime interval propagated by the z/j mass eigenstate. If 
this interval is large (s >> m~^) and time-like ( t > |x| ) then we can use the asymptotic expansion 

of the Bessel function 

/ 2 \i/2 



KAtms) - e^'""' (2.6) 

\7iimsJ 

to obtain from ( p.5|) the expression for the propagation amplitude at large time-like separation 

Kix, t- m,) = I , "^ ^ == I exp{-im^{tf-tY-{p,,-^,^)} (2.7) 



The phase factor of the amplitude ( p.7| ) is Lorentz-invariant and can be written in terms of the 
neutrino energy E and the time of flight t as measured from the lab frame as —irriiSi = —irrii (1 — 

^2)1/2 t = _,(^2 /^^) ^_ 

An extra factor of two appears on subtracting the phases at different spacetime points. A flavour 
eigenstates neutrino or kaon is observed at a single spacetime point and one should therefore compute 
the phase difference at the overlap of the two mass eigenstate wavepackets. This averaging over 
spacetime is done formally by considering the propagators of gaussian wavefunctions as opposed to 
plane waves , and the standard expression for the oscillation length is recovered. 

3 Conservation laws for long distance propagators 

A flavour eigenstate propagates as a linear combination of different mass eigenstates which are 

on-shell. The phase difference is obtained by some authors by assuming that the different mass 

eigenstates have common energy, and by some by assuming that they have a common momentum 

T6| . In this section we show that the linear combination of different mass eigenstates have neither 



the same energy nor the same momentum. We show that particles propagating over large distances 
(X >> P/m?) are constrained by the conservation laws at the vertex to be on shell. The on-shell 
condition is all that is needed to fix the phase difference and the oscillation length. 
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Consider a diagram with a propagator between vertices at spacetime points (xi,X2) with a 
number of external legs at the vertices. The amplitude is proportional to 

d'^Xi d!^X2 exp{—i'^qi-Xi) K{xi,X2) exp{i'^qf ■ X2) (3.1) 

■' i f 

Where qi are the incoming four momenta from the external legs at Xi and qj are the incoming four 
mommenta from the external legs at X2- Substituting the propagator 

d'p- ^— - (3-2) 

in (|3.1| ) and integrating over xi and X2 gives the usual energy mommentum conservation 5'^iJZi (li~p) 
and S^iJ^flf — p) at the vertices. When the proper distance between the two vertices s is larger 
than m~^, then the propagator used in (|3.1|) is the asymptotic form (|2.7| ), which in mommentum 



space can be written as |17 



K{x,,X2;s»m ') = ^^^^^^e" 



■ims 



f-=L= exp {-z(VP2 + m^){tf - U) + zP ■ (xf - x;)} (3.3) 



Substituting ( p.3|) in (|3.1| ) and integrating over xi we obtain 



f d^xi exp {-i{J2 Ei - VP2 + m2 )ti + i{J2 qi - P) ■ xi} 



= S{J2E,-VP^ + m^) <53(^q._p) (3.4) 

i i 

This implies that at a vertex energy and mommentum are conserved and in addition those particles 
which propagate without freely over distances larger than (|P|/r7i^) obey the mass shell constraint 
, E = v/p2 + m?. 

4 Path Integral method 

From the derivation of the propagation amplitude (|2.7| ) it is clear that there is no difference in the 
result on whether the propagating fields are bosons or fermions. The same result for the propagation 
amplitude is therefore also expected in relativistic quantum mechanics . In relativistic quantum 
mechanics the propagation amplitude K{mi, t) is computed from the classical action for a relativistic 
particle by the path integral method [l^. The amplitude for a particle of mass rrii to propagate 
from the spacetime point Xj to another spacetime point Xf can be written formally as, 

T^/ \ /"^r / M , im r^f ( dx'^ dx" ,\ , . ,, ^. 

K{xf,Xi;m) = J V[x{t)\ exp { — —J \Vf^u-T^-^ + ljdT} (4.1) 

For infinitesimal propagation between points Xj and Xj+i the quantum mechanical transition 
amplitude may be written as 

K{x,^,, X,) = C exp{-'^ i^^^%^ + 1) ^-} (4-2) 



The condition K{xj^i,Xj) = S'^i^Xj+i — Xj) as Ar -^ fixes the normahsation constant C = 
{m/2TTiAT) . In what follows we will not write out the explicit form for the normalisation factors 
as the final normalisation can be fixed from the boundary conditions. 

To carry out the path integral we divide each path into N —>■ oo infinitesimal segments of length 
Ar —>■ such that NAt = t the total path length. The amplitude for propagation from Xi to Xf 
is expressed as a product of the infinitesimal amplitudes K{xj,Xj^i) over the small straight line 
segments Ar, 

K{xf,Xi)= I dxi- ■ -dxN-i K{xf,XN-i)---K{xi,Xi) (4.3) 

By substituting the infinitesimal amplitudes ([4.2| ) in the expression ([4.3|) we have the amplitude 
K{xf,Xi) as a product of ordinary integrals of the infinitesimal amplitudes K{xj,Xj-i) over the 
intermediate spacetime points (xi, • • -xn-i)- 

K{xf,x,)=C' Jdxi---dxN-i exp{^if2(^J-^J+if + ^^A} (4-4) 



The Gaussian integrals in ( |4.4| ) can be performed sequentially and the result of the N—1 integrations 



IS 

' N 



j dxi---dxN-i exp {^^ I ^ (xj - Xj+i)M } = C^ ^^P^\2A^)\N^^'~^f'>J^ '^^'^^ 



Using ( [4.5|) in ( [4.4|) and writing A^Ar = r in the limit N —> oo and Ar ^ , we have the expression 
for the propagation amplitude K{xi, Xf) as a function of the proper time of r along some path, 

K(..,.„r) = {£-)' e.p{^ (i^i^ + r)} (4.6) 

where the normalisation constant has been fixed by the requirement that K{xi, Xf, r) = (5^(xj — Xf) 
as r — > 0. 

In relativistic quantum mechanics the proper time along a path r is an extra parameter that 
must be integrated over to include all paths including those which go backwards in the coordinate 
time t. The relativistic propagation amplitude from (|4.6| ) is 



r°° / m \'^ —im f (xi — XfY \ , ^ „, 

^(-/.-<) = /„ dr(—) exp— (L-^ + .j (4.7) 

The integral over r can be performed by using the formula ( |2.4| )to give the amplitude ( |4.7| ) as a 

Bessel function 

/ m \ 2 
K{xf,Xi) = ^—j - K^{ims) (4.8) 

which is of the same form as the field theory propagator ( p . 5|) its asymptotic form 

K(xf,Xi;s»m~^)c^i ) expi—ims) (4.9) 

\27rs/ 

is the same as (|2.7|). 

The expression given in ( |4.8| ) is the correct form of the propagator for particles which in general 
are not on-shell. To compute the propagator for strictly on-shell (real) particles, we and restrict 
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ourselves only to those paths which go forward in time and propagate over a fixed total proper 
time r = s = J{xf — Xj)^. The on-shell propagator is then obtained by inserting (2nT/mY^'^ S{t — 
Xf — Xi)"^) (the factor of (2nT/mY^'^ is for normalisation) in ([4 .71 ) to give, 

K{xf,Xi)on-sheU^ {Y~) (^xp{-ims) (4.10) 



Comparing ( [4.10| ) and ( [4. 91) we see that in the asymptotic limit the off-shell propagator becomes 



identical to the on-shell propagator. In other words 'virtual' particles become 'real' on large distance 
propagation. 

In non-relativistic quantum mechanics the proper time t = {tf — ti) the coordinate time inter- 
val. The non-relativistic propagator (which off-course is always on-shell) is obtained by inserting 
(27rr/m)^/^ 5{t — [tj — tif) in ( [4.7|) . The non-relativistic propagator thus obtained is the familiar 
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A™(x.,i.;x„i/) = (^-(^)"' e.p{-iMt,-U) (i - IpfJ)} (4.11) 



The amplitude ( 4.11 ) has the same form as the non-relativistic limit of the asymptotic propagator 
( p.7|) . As discussed earlier the phase factor exp — imt implies that in interference between two 
different mass states of non-relativistic particles (like Ki and Ks) the oscillations are periodic in 

{Am)t. 

5 Wavepacket analysis 

The propagator K{xf, Xi) given in (p?7| ) is the amplitude for the propagation of a particle localised 
at Xi (a delta function initial wave function) to be detected at Xf. The propagator of some general 
initial wavepacket \l/j„(x — Xj) is obtained by using the expressions for the delta function propagator 
given in ( p.7| ) and the superposition principle, 



K{xf — Xi) = I d X K{xf, x) "^inix — Xi) (5.1) 

If the initial wave function is a gaussian , 



2 (+ + \2 



fx-x-)^ it-t 
^in{x -Xi)=N exp{-t{Ea{t - U) + zPa • (x - Xi)} exp {- ^ ^ '' - , ^' } (5.2) 

where ax and at are the uncertainties is the initial position and time of production of the particle 
respectively and A^ is the normalisation constant. In the earlier wave-packet analyses of oscillation 
problem [^, the uncertainty in the time of production was neglected. In a covariant treatment 
both should be included. In the next section we show that the time uncertainty gives rise to a 
novel phenomenon of conversion without oscillations in many experimentally relevant situations. 
Substituting the expression for K{xf,x) given in ( |2.7D and the gaussian initial wavepacket ( |5.2D in 



JT]) and evaluating the integral by the stationary phase method [0 we obtain the expression for 
the propagation amplitude of a gaussian wavepacket , 

^(X, T; mj = (-^)^/^( .J", , )-'/'exp{-z^T + .P. ■ (X - v.T)} exp{-g^^^} 
47r|X|2 al + v^af Ea A[al + vlat) 

(5.3) 



where X = Xf — Xi and T = tf — ti are the space and time intervals propagated by the center of 
the gaussian wavepacket and Va = P^/Ea . The normahsation constant has been fixed such that 
/ dT d^X k{X,T;ma) k^{X,T;mb) = 1. The factor of (l/47r|XnV2 [^ (^ gho^g that the 
particle fiux decreases inversely with the sqaure of distance. Neutrino disappearance experiments 
look for evidence of depletion of a certain neutrino species over and above the expected inverse 
square decrease in fiux. In the following we will not display this factor (l/47r|Xp) in the probability 
expressions. The amplitude for the oscillation of one type of neutrino fiavour to another is obtained 
by substituting the propagation amplitudes for mass eigenstates (|5.3|) in , 

3 

A{a^P;X,T) = J2U^akima;X,T) Ul^ (5.4) 

a 

where a, /3 denote the fiavor eigenstates (z/e, t'^, Vr or i^°, K^) and the summation index a denotes 
a mass eigenstates (i^i,^'2,t'3 or Kl.Ks) ■ Uaa =< tt|a > and U^^ =< a\(3 > are the elements of 
the mixing matrix which relate the fiavor eigenstates with the mass eigenstates. The probabilty of 
fiavour oscilation as a function of space-time is the modulus squared of the amplitude (|L 



* |2 



P{a^P;X,T) = \Y,U^aK{ma;X,T) f/;^, 

a 

= T.\Upa\' \U^a\' + Y.Ui3aU*^aU;bUab k{X,T;ma) K^{X,T;mh) (5.5) 

In interference experiments the time of fiight of the particle is not measured and only the 
distance between the source and the detector is accurately known |T^ . The probability for the 
fiavour conversion as a function of distance is given by the time integral of (|5.5| ) , 

P{a^l3-X) = j dT\J2Uaak{ma-X,T) [/:^|2 

= Y.Ppa?\Uo.a? + Y^U^aKaUkU^, f dTK{X,T;ma) K^{X,T;m,) (5.6) 

The interference term given by the time-overlap of the propagation amplitudes (wave-functions) of 
different mass eigenstates can be evaluated for the gaussian propagator (|5.3| ) and is given by, 



Re JdT k{X,T-ma) k\X,T-m,) = - cos{{P^-P^) ■ X - {E, - E,y^^^^-^ ■ X} 



Xexp i [J^a J^b) ^o 8r^2 , „2^2^„2i ^^- 



where v = J{v'^ + vl)/2. In terms of the average momentum P and energy E and their respective 

differences APand AE the interference term ( p. 7] ) ,to the leading order in (AP/P) and (AE/E), 

reduces to the form, 

2 X 

- cos{ — {EAE - P ■ AP)} exp{-A} (5.8) 

where the exponetial damping factor, 

A = [AEf 2o- (|!) + {^f g (5.9) 



In general neither AE nor AP is zero and they depend upon how the state is prepared. For 
example if the mommentuni of the initial and the associated final state is measured then AP = 0. 
In the last section we have shown that for long distance propagation (s > m~^) P and E are not 
independent and are related by the mass shell constraint. The particular combination that appears 
in the phase difference turns out to be independent of the preparation and is fixed by the condition 
that each of the mass eigenstates be on shell, 

Ei'^ — P^ = m^ i = a,b 

^ 2EAE-2P-AP = Am^ (5.10) 

Substituting ( p.lO| ) in ( |5.8| ) we see that the interference term of ( |5.7| ) is given by , 

-cosi^j^X) e-^ , (5.11) 

The probability for the flavour conversion as a function of distance (|5.6|) is therefore, 
P{a^f3;X) = E - \Upa\'\U^a\^ 



1 ,,. ,..,. ,..,. ,. , Ami — rni^ 



+ E - \U,a K^U;, U^,\cos C \p '' X + 5) e-^ (5.12) 



where 5 = arg{Upa U^aU^b Uab)- 

We see that the standard oscillation formula which was obtained for relativistic particles [jl[ 
is actually valid at all energies. In other words although the formula ( 5.12|) is usually derived by 



taking the leading term in a {m? / P"^) series, our covariant calculation shows that there are actually 
no 0{m? /P'^) corrections to (|5.12|) . 

In the non-relativistic regime where P = {rria + mb)v/2, the standard kaon oscillation formula 
cos^nia — mb)T is recovered. 

There is no extra factor of two in the relativistic kaon oscillation formula as would have appeared 
without the wavepacket averaging . We see that although the two mass eigenstates have different 
proper times the wave packet overlap results in the average proper time appearing in the interference 
term. Therefore the interfernce phase is actually {ma — m,b)s = {ma — m,b){m,a + m,b)T/{Ea + Eb) = 
{rn^ — ml)T/2E. The average time prescription of P] can therefore be justified using the covariant 
propagator method. Other methods of showing that an extra factor of two does not appear in the 
flavour oscillations formula have been discussed in [HI. 



6 No muon oscillation 

The energy mommentum relation was employed in the last section to express the oscillation lenght 
in terms of the average mommentum of the mass eigenstates. We can use the conservation laws 
to give an expicit expression for the phase difference of the oscillation term. For example if the 
neutrinos are produced by pion decay vr^ -^ A*^^^, one can use the conservation laws to write the 
oscillation lenght of the neutrino in terms of m.,^ and rrinu and explicitly the validity of the formula 
for the phase difference (|5.11|) . In the pion rest frame the energy and the mommentum of the 
neutrino and the muon are, 

E. - 



ml — mf^ + ml 
2m^ 



E, 



ml + ml — ml 



2m 



■K 



IP I _ ((m^ - (m^ + m^f){ml - [m^ - m^)^))^/^ 



Using ( |6.1| ) we can check explicitly the expression for the neutrino phase difference, 



Ami 
2m.^ 



Am 
2m,T,F 

X Am^ 

A0, = (E,Ai^,-P,-AP,)— = -^X (6.2) 

Another application of the conservation laws is the question of associated particle oscillation 
raised by SWS [Q. For instance in the process vr"*" — >• //"''z/^, SWS state that since the muon 
mommenta depend on m,y, an interfernce of muon eigenstates corresponding to different Pij^{m,^) 
will result in the occurance of oscillation in space of the muon probability. Similarly in the reaction 
p p -^ A K , SWS claim that the associated oscillations of A probability in space ought to be seen. 
This claim has been disproved by Lowe et al and Okun et al, by differnt methods. We will see that 
the phase difference of m^ corresponding to the different z/ mass states actually vanishes . Using 
Ah we see that 



AE^ 



Ami 
2m.^ 



^ 2m^P ^ 

A(j)^ = (E^AE^-P^-AP^)^ = 0. (6.3) 

Okun et al state that muon oscillations do not take place, if the neutrino is not observed ' owing 
to the orthogonality of the neutrino eigenstates. That would imply that if there is a mixing of the 
neutrinos with some heavy states , which would make the low-lying mass neutrinos non-orthogonal, 
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then muon oscillations will take place. This would have been powerful method for testing for the 
orthogonality of the observed neutrinos. However from ( |6.3| ) we see that the form of the phase 
difference is such that the conservation laws ensure that whether the v is measured or not there 
is no EPR correlation to the [i phase difference. There is an EPR type correlation of E^;^_ and P^ 
with the neutrino mass eigenstates but the phase difference vanishes . Therefore there will be no 
muon oscillations In the linear combination of different mass eigenstates, niether !S.E nor AP is 
zero. Assuming that either of them is zero would lead to incorrect results. 

7 Conversion without oscillation 

The suppression factor A which goes with the oscillation term has some interesting new implications. 
If the parameters of the experiment are such that A becomes large then no oscillations is spacetime 
can be observed. What can be observed is a constant conversion probability. 

For example, the survival probability of a z/^ after propagating over a distance X in the large A 
limit is, ( for two neutrino flavours (z/g, v^ with the mixing matrix elements t/gi = f^^2 = cosQ and 
t/e2 = -Vy.i = sin6), 

P{v^ -^ u^; X) = sm^e + cos^O + 2sin'^ecos'^e cos{-^X) e"^ ~ sm^O + cos^O (7.1) 

and the expression for the conversion probability when A is large is, 

P{v^ -^ Ve, t) = 2sin^e cos^e (1 - cos{^—X) e"^) ^ 2sin^e cos^O (7.2) 

It has been noted ealier [Q one condition A must be small which implies that X must be smaller 
than the coherence length L^oh = (4A/2aP^/Am^) . We shall show below that the constraint 
Lose < Lcoh is not a sufficient condition to ensure the occurrence of flavour ocillations in space. In 
the analysis of refs |^ only the uncertainty in the initial position a^ was considered, in our analysis 
we have included the contribution of at the uncertainty in time at which the neutrino was produced 
, which as it turns out, makes the larger contribution to the suppression term A. This happens when 
the neutrinos are produced from long lived resonances. The neutrino wavepacket which is produced 
has a spread in time with width at which cannot be smaller than the lifetime of the resonance whose 
decay produces it. The dominant contribution in that case arises from the first term of A, 

where we have equated the initial time uncertainity at with r - the lifetime of the resonance that 
produces the neutrino, and we have ignored the spatial spread of the wave-packet ax which in 
most experiments is many orders of magnitude smaller than at- In terms of the oscillation length 
Lose = 4:7iE/{Am'^) the suppression factor A. In order to observe oscillations the source-detector 
distance L must be larger than Lose- Which that the minimum source-detector distance in order 
that neutrino oscillations be observed is given by condition, L^m = (vr)-\/2T. We shall show below 
that most neutrino experiments do not satisfy this criterion for observability of space oscillations. 
On the other hand spatial oscillations of fiavour can be observed in the $ and B factories owing to 
the large width of and T resonances. 
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For relativistic particles produced from long lived resonances the formula for the conversion 
probability which should be fitted with the experiments is , 

P{^a ^ ^/3) = 2 "'^2^ y- - ^^"( E ^ ''''^^'^ E ^^ ) ^ ^^ 

where Am^ is the mass square difference in eV"^ ,L is the detector distance in m (km), r is the 
lifetime of the parent particle in the lab frame in m {km) and E is the energy in MeV (GeV). 

8 Constraints on Am^ and sin^26 from experiments 

We shall now fit the oscillation formula (774) with the results of some terrestrial neutrino oscillation 



experiments , classified according to the type of the source of neutrinos. 

Stopped Muons : The LSND experiment |^ has searched for i?^ ^ z/g in a beam of anti-muon 
neutrinos produced from stopped muons. The minimum uncertainty in the time of u^^ production is 
the muon lifetime r^ = 658.65m. The average neutrino energy E^ ~ 30MeV and the detector is at 
a distance L = 30m from the source. In order that neutrino oscillations from muon decay neutrinos 
be seen the minimum source detector distance must be Lmin = V^t^t^ = 2927.6m. This means that 
the conversion probability does not oscillate as a function of distance. Fitting the reported range 
of conversion probability P{iy^ —>■ Ve) = (0.31 + 0.11 — 0.1 ±0.05) x 10"^ with the formula (ff^) , the 
regions of allowed Am^ and sin'^26 which gives rise to this probability band are shown as dotted 
band in Fig. 1. The same range of conversion probability when fitted with the standard oscillation 
formula gives the shaded region shown in Fig.l. We have averaged the conversion probability by 
assuming a gaussian distribution of neutrino energy, with average energy SOMeV^ and a half width 
of lOMeV. The asymptotic bounds are as follows. If Am^ > ( 2^E/t^) = 2.54 x lO-^el/^ then 
mixing angle constrained in the range 0.3 x 10~^ < si'n'^26 < 0.9 x 10^^. For large mixing angles 
sin^2^ ~ 1 the mass square difference lies in the range, 2 x lO'^el^^ < Am^ < 3 x lO^^el^^. This 
bound is two orders lower than what is obtained by the use the standard oscillation formula 0. 

Stopped Pions : The Karmen experiment |jll[ searches for u^ —>■ u^, conversion in a beam of 
u^ produced by the decay of stopped pions . The time at which the z/^ is prduced cannot be 
localised to less than the pion hfetime Tt,+ = 7.804m. The u^ energy is E,^ = 29.8Me\^. In order 
for oscillations to be observable the detector distance must be larger than L = \pl'nTT^^ = 35m. 
Since the detector distance is 17.5m, it is not possible to observe spatial oscillations in the Karmen 
experiment. The experimental bound -P(t'^ — >■ z/g) < 3 x 10~^, fitted with the conversion probabilty 
formula (O) constrains sin'^2e < 0.6x10-^ if Am^ > {2V2E/t^+) = 2 eV^. In the sin^2e ~ 1 limit 



Am"^ < 0.16ey^. We have averaged the conversion probability by assuming a gaussian distribution 
of neutrino energy, with average energy SOMeV and a half width of lOMeV. The regions of the 
parameter space allowed by the covariant wavepacket formula and by the standard formula are 
shown in Fig. 2. Since the pion decay length is of the same order as the experimental baseline, the 
improvement in sensitivity to Am^ is marginal compared to LSND. 

High energy Pions or Kaons: When high energy neutrino beams are obtained from the decay 
of a pion or kaon in flight the both r and E are larger by the Lorentz factor 7 and consequently A 
remains the same as that with stopped pions or kaons. In the BNL E776 experiment |T^ most of the 



neutrinos are produced in vr^ decay and have an average energy E ~ ^Gev and the detector is at a 
distance L = 1km . In the lab frame the neutrino time uncertainty is r = r^+ (i?^+/m^+). Taking the 
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mean tt^ energy to be lOGeV , the lifetime of the pions in the lab frame is r = 7.804(10/0. 135)?7i = 
0.578km. Fitting the E776 experimental limit -P(t'^ — > i^e) < 1.5 x 10^'^ with the conversion 
probability formula (^, we see from Fig.3 that sin^29 < O.SxlO'^ if Am^ > {2V2E/t^+) = 5 eV"^. 
In the sin'^2e ~ 1 hmit , Am^ < 0.2e\/2 . 

In the CCFR |]12| experiments most of the neutrinos are produced in K^ decay and have an 



average energy E ~ lAOGev and the detector is at a distance L = l.Akm . In the lab frame the 
neutrino time uncertainty is r = tk+{Ek+ /TnK+). Taking the mean K^ energy to be 600GeV^ we 
see that the lifetime of the kaons in the lab frame is r = 3. 7(600/0. 493)m = A.blAkm. Fitting the 
CCFR experimental limit ^(t'^ -^ Ve) < 0.9 x 10^^ with the conversion probability formula ( |7.4| ) 
we see from Fig.3 that sin'^2e < 1.8 x lO'^ if Am^ > {2V2E/tk+) = 76 eV^. In the sin'^2e ~ 1 
limit , Am^ < 1.2ey^ 

Reactor neutrinos: Reactor neutrinos are produced primarily from the beta decay of ^35^/ 
and "^"^^Pu which have a half life ~ lO^sec |T3|. The minimum detector distance for observ- 
ing neutrino oscillations is lO^^m which means that it is impossible to detect spatial oscillations 
with reactor neutrinos. The upper bounds on conversion probability in the Bugey experiment is 
P(Pe -^ ^pL or Ut) < 0.075. The largest detector distance L = 95m and the average neutrino energy 
is E' ~ 5MeV. The conversion probability formula is independent of Am^ if Am^ > lO^^eV^^. Fig 3. 
shows that in this limit the reactor experiments rule out sin'^29 > 0.15. In the standard oscillation 
formula the oscillatory term averages to zero when Am"^ > {4ttE/L) . Using the wavepacket formula 
( [7.41 ) however we see that the oscillatory term is is exponentially damped at much lower values of 
Am^ . If neutrino sources have a large decay time which is the case in most experiments , the 
conversion probabilty is sensitive to mass differences Am? > {2\/2E/t). For this reason the LSND 
experiment which uses neutrinos from muon decay has two orders of magnitude more sensitive than 
KARMEN, BNL-E776 etc where the neutrinos are from vr and K decay. This is evident from the 
combined plot of these experimental results, shown in Fig.3 fitted with the covariant wavepacket 
formula (|7.4|) . The reactor experiments where r ~ lO^s are in fact sensitive to values of Am"^ as 
low as 10~^eV'^. But the probability measurement in reactor experiments is poor compared to the 
accelerator experiments which is why they rule out only a small region of parameter space. 

Meson oscillations at B and $ factories: The mesons produced by or T(4S') decay are non- 
relativistic and the expression for the damping factor is A = 2 (Am r)^. For the kaons produced 
from from decays, Am^ = 3.5 x 10~^'^MeV and the lifetime is r^ = 0.225MeV~^, and the 
damping factor A = 2.4 x 10"^'* is negligible. For the case of BB oscillations produced from T(45') 
decays, the Am^ = 3.53 x 10"^° MeV and the lifetime of T(4S') is tt = 0.04:2MeV-^ and again the 
damping factor, A = 8.6 x 10"^^ is negligible. This is important for tests of EPR [J19l and Quantum 



Mechanics [20| where the measurement of the oscillatory term of the probability is essential. 



9 Conclusions 

The covariant formulation gives two new results compared to the standard treatment. We use the 
condition that the spread of the neutrino wave function along the time axis is given by the lifetime 
of the particle whose decay produces the neutrino. When neutrinos are produced from long lived 
particles like muons , the interference term in the conversion probability formula vanishes for much 
smaller values of the neutrino mass square difference compared with the standard formula. The 
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reason this happens can be understood with the following picture. When the spread of the neutrio 
is large along the time axis, its spread along the energy axis is small. The neutrino wavefunctions 
for different masses are then energy eigenstates which are mutually orthogonal. The interference 
term which is the overlap of the wavefunctions of diffferent mass states therefore vanishes when 
these states become orthogonal. 

The second new result is that the for the secondary particles, like muons in the decay n ^ fiu 
there is no spacetime variation of probability independent of whether the neutrinos are measured 
or not. Although the muon wavefunction is a linear combination of distinct energy and momentum 
states corressponding to the different u mass states, the phase difference for the secondary par- 
ticles vanishes. Therefore there is no EPR type oscillation of the secondary particles probability 
distribution. This is true independent of whether the primary oscillations of the z/'s is measured or 
not. This conclusion is results from the application of the correct conservation laws according to 
which the different mass states of u have neither energy nor momentum in common as is sometime 
assumed in the standard derivations. 

Acknowledgments I thank Terry Goldman, Walter Grimus and Marek Nowakowski for useful 
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Figure 1: Lsnd n decay at rest experiment ||^ allowed regions with the standard formula (between 
the dotted curves)and the covariant oscillation formula (continuous curves). 



15 




-2.5 



-1.5 -1 

Log(sin'^2 2 theta) 



-0.5 



Figure 2: Karmen vr decay at rest experiment ||TT| allowed regions with the standard formula 
(dotted line) and the covariant oscillation formulas (dashed line). 
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Figure 3: The region between the continuous hues is allowed by the LSND /j, decay in flight experi- 
ment [§ .The region between the dotted lines is allowed by the LSND vr decay in flight experiment 



21| . Region ruled out by E776 |]T0[ is above the dashed-dotted curve and by Karmen |[TT[] is above 
dashed curve. Region above the topmost dashed curve is ruled out by CCFR UTT 
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